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INTRODUCTION 

Some  applications  of  relativistic  electron  beams  (REB)  necessitate  that 
they  stably  propagate  some  distance  through  air.  Whether  the  beam  will  propagate 
or  not  depends  on  the  pressure  of  the  air  into  which  it  is  fired.  The  goal  of 
many  REB  experiments  is  to  determine  the  pressure  most  conducive  to  beam 
propagation. 

Low  pressure  environments,  into  which  an  REB  may  be  fired,  can  be  created 
by  various  means.  The  most  practical  means,  experimentally,  is  to  use  a  drift 
tube  whose  interior  pressure  is  adjustable.  The  beam  is  then  fired  into  the 
tube  and  monitored  as  it  passes  various  points  along  the  length  of  the  chamber. 

The  drift  tubes  used  in  propagation  experiments  generally  have  cylindrical 
symmetry,  and  are  sometimes  made  of  highly  conductive  metals.  Hence,  because  of 
the  conducting  boundary,  firing  an  REB  into  a  metallic  drift  tube  may  not  be  the 
same  as  firing  it  into  the  unbounded  atmosphere. 

In  this  investigation,  we  endeavor  to  determine  (by  solving  Maxwell's 
equations  with  appropriate  boundary  conditions)  a  minimum  drift  tube  diameter 
such  that  conducting  boundaries  have  a  negligible  affect  on  relativistic  electron 
beam — plasma  interactions. 
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ASSUMPTIONS  AND  EQUATIONS 


As  is  well  known,  Maxwell's  equations,  in  Gaussian  units,  for  material  media 

are: 

7*1-  4irp  (1) 

7*1-0  (2) 

7x1-  -i||  (3) 

'>*■7  I  +  7i!  (4> 

To  obtain  the  appropriate  wave  equation,  apply  the  7x  operator  to  equation  3 
and  employ  the  identity 

7  x  7  x  F  -  7  (7*F)  -  72  F 


to  get 


72  I  -  7  (7*1)  --£-7x1 


Substituting  into  equation  5,  the  expressions  for  7  •  E  and  7  x  B  yields 
<72  ?  1  32£  4tt  9J 

7  E  -  o.  JO  "  T  (6) 

To  investigate  the  problem,  we  need  solve  only  the  z  component  of  equation  6. 

To  do  so,  we  make  the  assumption  of  cylindrical  symmetry.  That  is,  we  assume 

that  Eq,  B^,  and  Bz  are  zero,  and  that  the  remaining  components  of  E  and  B  are 

independent  of  the  azimuthal  coordinate  In  addition,  we  assume  that  Ohm's  law, 
*►  *► 

J  -  cjE,  is  valid  in  Che  plasma.  With  the  above  assumptions,  the  z  component  of 
equation  6  becomes 

3nF--L32F 

r  ^r  3r  z  3z2  z  c*  Tt  z  c2  St2  z 


where  J'  represents  the  REB  current  density, 
z 
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We  simplify  equation  7  via  the  transformation 


Now 


_3_  _3_ 

3t  ^  3t 


(8a) 

(8b) 


(9a) 


and 

_3_  ^  _3_  1  _3_ 

3z  3Z  “  c  3t 


(9b) 


We  now  make  the  so-called  frozen  beam  approximation;  i.e.,  we  assume  that  the 
fields  translate  rigidly  with  the  beam.  For  any  given  field  component  F,  the 
frozen  beam  approximation  is  equivalent  to  the  mathematical  statement 


DF  3F  3F 

Dt  "  3t  3z 


0 


or 


3F  3F 

3t  “  "  v  3z 


(10) 


(ID 


where  v  is  the  velocity  of  the  beam  electrons .  The  energy  of  the  beam  electrons 
is  generally  greater  than  1  MeV.  Hence,  we  can  further  assume  that  v  *  c.  Thus, 


3E  3E 

_ z  z 

3t  c  3z 


(12) 


Substitution  of  the  operator  identities  9a  and  9b  into  equation  12  yields 


3F. 

3E 

3E 

z 

z 

z 

3t  “ 

3t 

c  3Z 

(13) 


Equation  13  is  valid  only  if 
3E 

c  -r#-  2  0  (14) 


That  is,  with  respect  to  the  coordinates  Z  and  T,  the  frozen  beam  approximation 
is  equivalent  to 


3E 

TT  '  0  (15) 
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In  view  of  equation  15,  the  identities  9a  and  9b  become 

1^1 

St  St 


In  terms  of  the  new  coordinates,  equation  7  can  be  written  as 


(16b) 


-  T-  r  E  -  J-  a  E 
r  Sr  Sr  z  c2  St  z 


4tt  Sp  47T  ^^2 
c  St  c2  St 


An  REB  propagating  through  air  generates  a  plasma,  resulting  in  space  charge 
neutralization  (SCN).  We  account  for  SCN  by  letting 


p  *  —  exp  (-  4ttot) 


Substituting  equation  18  into  17  yields 


13  3  _  4ir  S  47rS  , 

r  Sr  r  Sr  Ez  c2  St  z  c2  St  ^z 


where 


J'*  ■  J'  [l  -  exp  (-4ttot)] 
z  z 


We  assume  that  the  conductivity  of  the  drift  tube  is  large  enough  to  be  con¬ 
sidered  infinite.  Continuity  of  Ez  across  the  plasma — wall  interface  requires 


M 


where  Rw  represents  the  inner  radius  of  the  drift  tube. 

Since  the  beam  electrons  are  highly  relativistic,  Ez  is  approximately  zero 
at  the  head  of  the  pulse.  That  is, 


E_  (r,o)  -  0 

Z 
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The  0  component  of  equation  3,  in  the  new  coordinates,  is 

h  Ez  (r-T) " "  cir  Bs  (r»T) 

Evaluating  equation  23  at  r  *  0  yields 
3 


(23) 


37  EZ  (r’T) 


r-0 


(24) 


Equations  19,  21,  22  and  24  uniquely  determine  Ez  (r,T)  inside  the  drift 


tube  chamber. 
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BEAM  PROFILE  AND  PLASMA  CONDUCTIVITY 

The  beam  is  assumed  to  have  the  following  profile: 


J°  (T/Tji/d  +  r2/a2)2 

TiTi 

(25a) 

J°  /(I  +  r2/a2)2 

1  T  1  T2 

(25b) 

0 

T  >  T 

(25c) 

where  J°  and  a  represent  the  on-axis  beam  current'  density  and  Bennet  radius, 
respectively.  Also, 

j£  «  i£/7ra2  (26) 

where  1°  is  the  on  axis  beam  current. 

We  assume  that  the  plasma  conductivity  (O)  is  generated  with  a  radial  profile 
similar  to  that  of  the  beam.  That  is,  we  let 

a  *  a0/(l  +  x2 fa.2)2  (27) 

where  afl  is  the  on-axis  conductivity.  We  employed  existing  air  chemistry  models 
(Ref.  1)*  to  numerically  obtain  the  on-axis  conductivity.  Figure  1  is  a  graph  of 
on-axis  conductivity  versus  distance  (in  seconds)  back  into  the  pulse.  The  beam 
parameters  used  are  indicated  at  the  lower  right-hand  corner  of  the  graph.  The 
rise  length  of  the  pulse  is  AL(=  Tjc)  and  its  length  is  L(«  T2c)  • 


1  Dreyer.  ’  A.,  Electric  and  Magnetic  Fields  of  An  Intense  Pulse  of  Relativistic 

.’sot;  <s  Propagating  Through  Air,  AFIT/DS/PH/79-1,  AFIT,  Wright-Patterson  AFB 
..ao,  1979. 

*  Johnston,  R.,  private  communication. 
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Figure  1.  On-Axis  conductivity  versus  distance 
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NUMERICAL  ALGORITHM 


We  solve  equation  19  using  an  implicit  finite  difference  scheme  developed  by 
Crank  and  Nicholson  (Ref.  2).  Following  this  scheme,  equation  19  becomes 


£jF[K)ri-K)(Ht  va 


?n+l 

-zi 


+  i  —  6  r .  6  En . 
2  r  i  r  zi 


4tt  _L_  /  '*n+l  '*n  ^ 

"  cz  At  V  zi  "  Jzi  / 

The  boundary  conditions,  equations  21,  22  and  24, become 


(28) 


where 


:°.  -  o 

Zl 

(29) 

:n^  *  o 

zl 

(30) 

,n+l  =  En+1 

(31) 

zi+1  zi 

■  »  nAt 

(32) 

and 


r  *  iAr 


R  -  IAr 
w 


(33) 


(34) 


The  step  sizes.  At  and  Ar,  were  assigned  the  values  5.00  x  10-12  s  and  0.01  cm, 
respectively. 

The  z  component  of  the  plasma  current  density  ^Jz  j  is  calculated  using  Ohm's 
law.  That  is, 


jn+i  „  0n+l  £n+l 
zi  i  zi 


(35) 


2.  Carnahan,  B. ,  Luther,  H.  A.,  Wilkes,  J.  0.,  Applied  Numerioal  Methods , 
J.  Whiley  &  Sons,  Inc.,  New  York,  1969. 
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RESULTS  AND  CONCLUSIONS 


The  results  of  our  computer  runs  are  summarized  in  Table  1  and  Figure  2. 

The  beam  parameters,  with  which  equation  28  was  solved,  are  given  in  the  lower 
right-hand  corner  of  this  figure.  The  plasma  currents  in  this  table  and  figure 
were  evaluated  at  the  top  of  the  current  rise  (t*Tj)  for  three  different  pressures, 
given  in  units  of  Torr. 


TABLE  1.  RESULTS  OF  COMPUTER  RUNS 


P-760T 


P-76T 


P-7.6T 


R  (cm) 
w 

I 

(KA) 

I 

P 

|  (KA) 

I 

1  P 

!  (ka) 

3.0 

3.3 

5.0 

6.5 

8.0 

9.5 

10.0 

10.7 

12.0 

11.7 

15.0 

12.7 

20.0 

13.9 

35.0 

15.7 

50.0 

16.7 

3.8 

8.5 

7.1 

14.0 

10.0 

17.9 

11.2 

19.3 

12.1 

20.3 

13.0 

21.3 

14.1 

22.3 

15.9 

23.8 

16.8 

24.5 

The  curves  in  Figure  2  indicate  that,  at  each  pressure,  the  plasma  current 
increases  significantly  out  to  Rw  ■  20.0  cm.  Thereafter,  out  to  Rw  *  50.0  cm, 
the  plasma  current  increases  only  a  few  percentage  points.  That  is,  the  conducting 
boundaries  significantly  influence  the  plasma  current  when  the  drift  tube  radius 
is  less  than  10  Bennet  radii.  Thus,  for  a  beam  with  a  2  cm  Bennet  radius,  the 
minimum  drift  tube  radius  beyond  which  wall  effects  are  insignificant  is  10 
Bennet  radii. 

On  the  basis  of  our  results,  it  is  suggested  that  regardless  of  beam  size 
wall  affects  are  negligible  when  the  ratio  of  drift  tube  radius  to  Bennet  radius 
is  at  least  10.  That  is,  drift  tube  propagation  experiments  may  simulate  propaga¬ 
tion  of  beams  through  the  air,  provided  that  the  ratio  of  beam  radius  to  drift 
tube  radius  is  greater  than  10. 
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P  *7.8T 


P=7B.0T 


P  =  780.0T 


l„  =  30M 

T#  = 1 . 5MeV 
a  =2. Ocm 

A  L  *600cm(20ns) 
L*412Scm(137.5ns) 


S.O  tO.O  1S.0  20.0  25.0  30.0  35.0  40.0  45.0  50.0 

Rw(CM) 

Figure  2.  Results  of  computer  runs. 
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